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Abstract 

This  paper  reports  on  a  visual  study  of  the  CO2  bubble  behavior  in  the  anode  flow  field  of  an  in-house  fabricated  transparent  Direct 
Methanol  Fuel  Cell  (DMFC),  which  consisted  of  a  membrane  electrode  assembly  (MEA)  with  an  active  area  of  4.0  x  4.0  cm2,  two  bipolar 
plates  with  a  single  serpentine  channel,  and  a  transparent  enclosure.  The  study  reveals  that  at  low  current  densities,  small  discrete  bubbles 
appeared  in  the  anode  flow  field.  At  moderate  current  densities,  a  number  of  gas  slugs  formed,  in  addition  to  small  discrete  bubbles.  And  at 
high  current  densities,  the  flow  field  was  predominated  by  rather  long  gas  slugs.  The  experiments  also  indicate  that  the  cell  orientation  had 
a  significant  effect  on  the  cell  performance,  especially  at  low  methanol  flow  rates;  for  the  present  flow  field  design  the  best  cell  performance 
could  be  achieved  when  the  cell  was  orientated  vertically.  It  has  been  shown  that  higher  methanol  solution  flow  rates  reduced  the  average 
length  and  the  number  of  gas  slugs  in  the  flow  field,  but  led  to  an  increased  methanol  crossover.  In  particular,  the  effect  of  methanol  solution 
flow  rates  on  the  cell  performance  became  more  pronounced  at  low  temperatures.  The  effect  of  temperature  on  the  bubble  behavior  and  the 
cell  performance  was  also  examined.  Furthermore,  for  the  present  flow  field  consisting  of  a  single  serpentine  channel,  the  channel-blocking 
phenomenon  caused  by  CO2  gas  slugs  was  never  encountered  under  all  the  test  conditions  in  this  work. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

It  has  recently  been  shown  that  the  opportunity  in  terms 
of  market  size  and  growth  rate  for  small  fuel  cells  for 
portable  electronic  devices  is  commercially  compelling  as 
compared  with  batteries  [1].  Since  methanol  has  a  higher 
energy  density  under  ambient  conditions  and  can  be  pro¬ 
cessed  and  stored  more  easily  than  a  gas  fuel,  a  liquid  fed 
direct  methanol  fuel  cell  (DMFC)  is  becoming  a  potential 
power  source  potable  applications  [2,3].  For  this  reason,  the 
DMFC  is  now  being  studied  by  a  number  of  university  re¬ 
searchers  and  industry  engineers  all  over  the  world  [4-20]. 

The  DMFC  consists  of  a  proton  exchange  membrane  with 
an  anode  and  a  cathode  catalyst  layer  on  each  side.  Gas  dif¬ 
fusion  layers,  usually  made  of  carbon  paper  or  carbon  cloth, 
are  used  to  cover  the  catalyst  layers  and  form  the  membrane 
electrode  assembly  (MEA).  The  MEA  is  then  sandwiched 
between  two  current  collectors  that  have  some  flow  chan¬ 
nels  machined  into  the  surface  for  the  supply  of  fuel  and 
oxidant.  Typically,  the  DMFC  is  operated  at  a  temperature 
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lower  than  100  °C.  At  the  anode,  methanol  is  electrochem- 
ically  oxidized  to  CO2,  while  at  the  cathode,  oxygen  is  re¬ 
duced  to  water:  the  electrochemical  reactions  taking  place 
in  the  cell  are  described  as  follows: 

CH3OH  +  H2O  -+  CO2  +  6H+  +  6e_  (anode  reaction) 

(1) 

|C>2  +  6H+  +  6e_  -+  3H2O  (cathode reaction)  (2) 

which  can  be  combined  to  give  the  overall  reaction  as: 
CH3OH  +  f  02  -*  2H20  +  C02  (3) 

Thus,  the  overall  cell  reaction  is  the  electro-oxidation  of 
methanol  to  CO2  and  water.  Apparently,  all  the  reaction 
products,  gas  CO2  at  the  anode  and  liquid  water  at  the  cath¬ 
ode,  should  be  removed  from  the  electrode  structure  and 
cell  as  efficiently  as  possible  to  maintain  an  effective  con¬ 
tinuous  reaction.  Like  the  water  management  at  the  cath¬ 
ode,  the  efficient  removal  of  CO2  at  the  anode  is  one  of 
the  most  important  research  issues  in  the  development  of 
DMFCs. 

In  the  DMFC,  methanol  solution  is  fed  into  the  anode 
flow  field  and  diffuses  to  the  catalyst  sites  through  the  gas 
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diffusion  layer,  while  the  reaction-produced  gas  CO2  at  the 
catalyst  layer  transport  backward  into  the  anode  channels 
through  the  gas  diffusion  layer.  As  such,  a  counter  liquid-gas 
two-phase  flow  takes  place  in  the  diffusion  layer.  Since,  the 
quantity  of  CO2  increases  with  increasing  the  current  den¬ 
sity,  at  high  current  densities,  the  volume  fraction  of  CO2 
in  the  anode  flow  field  becomes  rather  large.  Under  such  a 
situation,  CO2  removal  from  the  catalyst  sites  is  critical  to 
ensure  availability  of  an  adequate  surface  area  for  methanol 
oxidation.  To  this  end,  CO2  has  to  be  removed  effectively 
from  both  the  diffusion  layer  and  the  anode  flow  field.  Oth¬ 
erwise,  CO2  gas  slugs  accumulated  in  the  anode  flow  field 
may  not  only  hinder  the  methanol  transport  to  the  catalyst 
sites,  but  also  cause  a  higher  pressure  in  the  anode  flow  field, 
leading  to  an  increase  in  methanol  crossover.  Thus,  research 
on  two-phase  flow  characteristics  of  gas  CO2  and  methanol 
solution  in  the  anode  flow  field  is  essential  for  the  design 
and  optimization  of  a  DMFC. 

Relatively,  few  papers  have  been  reported  on  the  study 
of  CO2  bubble  behavior  in  the  anode  flow  field  of  DMFCs. 
Mench  et  al.  [21]  examined  gas  bubble  growth  and  ejection 
from  the  backing  layer/flow  channel  interface  region  with 
video  microscopy  and  observed  discrete  bubbles  of  the  or¬ 
der  of  0. 1-0.5  mm  evolving  from  various  locations  within 
the  backing  layer.  Argyropoulos  et  al.  [22]  studied  the  CO2 
flow  characteristics  in  a  DMFC  with  the  aid  of  a  high-speed 
video  camera  and  found  three  flow  patterns  including  bub¬ 
bly,  slug  and  annular  flow  in  practical  fuel  cell  operation 
conditions.  In  their  experiments,  they  tested  two  different 
flow  field  designs:  one  consisted  of  parallel  channels  hav¬ 
ing  a  cross-section  area  of  2.0  x  2.0  mm2  and  a  length  of 
30.0  mm;  the  other  consisted  of  cross  flow  channels  having 
the  same  size,  but  including  a  triangular  enlarging  inlet  and 
outlet  sections.  They  studied  the  effects  of  the  current  den¬ 
sity  and  the  methanol  solution  flow  rate  and  found  that  gas 
slugs  at  low  flow  rates  and  high  current  densities  blocked 
the  flow  channels.  Scott  et  al.  [23]  visually  investigated  the 
CO2  gas  evolution  and  flow  behavior  with  flow  beds  based 
on  stainless  steel  mesh.  In  their  works,  a  number  of  the 
flow  designs,  based  on  stainless  steel  mesh,  showed  promis¬ 
ing  behavior  in  terms  of  gas  removal  characteristics  and 
electrical  performance.  Bewer  et  al.  [24]  reported  a  novel 
method  simulating  two-phase  flow  in  a  DMFC  using  an 
aqueous  H2O2  solution.  The  influence  of  the  flow-field  on 
the  bubble  formation  and  on  the  flow  homogeneity,  as  well 
as  the  influence  of  manifold  on  the  flow  homogeneity  was 
investigated. 

In  this  paper,  we  present  a  visual  investigation  of  CO2  gas 
bubble  behavior  in  a  transparent  DMFC  with  its  anode  flow 
field  consisting  of  a  single  serpentine  channel.  The  in  situ 
taken  images  show  that  there  exist  various  CO2  gas  bubble 
flow  patterns  at  different  current  densities.  Effects  of  various 
operation  parameters,  including  methanol  solution  flow  rate, 
cell  orientation,  and  cell  operating  temperature,  on  the  cell 
performance  and  the  two-phase  flow  characteristics  in  the 
anode  flow  field  were  examined. 


2.  Experimental 

2.7.  Transparent  cell 

A  transparent  DMFC  was  designed  and  fabricated  for 
this  visualization  study,  as  schematically  shown  in  Fig.  1. 
The  MEA  detailed  in  the  subsequent  paragraph,  was  sand¬ 
wiched  between  two  bipolar  plates  with  a  gasket  onto  either 
side  of  the  MEA.  This  assembly,  including  the  bipolar 
plates  and  MEA,  was  clamped  between  two  enclosure  plates 
by  eight  M8  screw  joints,  each  having  a  torque  of  about 
3  Nm. 

The  MEA,  fabricated  in  this  work,  had  an  active  area  of 
4.0  x  4.0  cm2  and  consisted  of  two  single-side  ELAT  elec¬ 
trodes  from  E-TEK  and  a  Nation®  membrane  115.  Both 
anode  and  cathode  electrodes  used  carbon  cloth  (E-TEK, 
Type  A)  as  the  backing  support  layer  with  30%  PTFE 
wet-proofing  treatment.  The  catalyst  loading  on  the  anode 
side  was  4.0  mg  cm-2  with  unsupported  [Pt:Ru]  Ox  (1:1 
a/o),  while  the  catalyst  loading  on  the  cathode  side  was 
2.0  mg  cm-2  using  40%  Pt  on  Vulcan  XC-72.  Furthermore, 
0.8  mg  cm-2  Nation®  was  applied  onto  the  surface  of  each 
electrode.  The  Nation®  membrane  115,  before  being  used, 
was  first  cleaned  following  a  standard  procedure:  (i)  boil¬ 
ing  membrane  in  5  wt.%  H2O2  solution  of  80  °C  for  1  h; 
(ii)  rinsing  with  DI  water  of  80  °C  for  1  h,  (iii)  boiling 
membrane  in  0.5  M  H2SO4  solution  of  80  °C  for  1  h,  and 
(iv)  rinsing  with  DI  water  of  80  °C  for  1  h.  Finally,  the 
MEA  was  formed  by  hot  pressing  at  135  °C  and  5  MPa 
for  3  min. 

To  avoid  corrosion,  the  bipolar  plates,  shown  in  Fig.  2, 
were  made  of  316  stainless  steel  [25-28]  plates  with  a  thick¬ 
ness  of  2.0  mm.  As  can  be  seen  in  Fig.  2,  the  rectangular 
bipolar  plate  consisted  of  two  portions,  the  channel  area  and 
the  extension  area.  The  channel  area  acted  as  the  distributor 
for  supplying  fuel  and  oxidant  to  the  MEA,  in  which  a  sin¬ 
gle  serpentine  channel,  2.0  x  2.0  mm2,  was  machined  by  the 


MEA 

'Gasket 
Bipolar  plate 

Transparent  enclosure 


Fig.  1.  Transparent  DMFC. 
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Fig.  2.  Flow  field  and  the  current  collector. 


wire-cut  technology.  The  width  of  the  ribs  was  1.5  mm.  The 
serpentine  channel,  having  a  total  length  of  420  mm,  con¬ 
sisted  of  11  serially-linked  horizontal  segments  (35.0  mm 
long)  and  10  vertical  segments  (3.5  mm  long).  The  exten¬ 
sion  area  of  the  bipolar  plates  served  as  a  current  collector. 
In  addition,  a  tape  heater  was  attached  to  the  extension  area 
to  adjust  the  cell  operating  temperature  to  a  desired  value 
during  the  experiments. 

For  the  purpose  of  visualization,  the  two  enclosure  plates, 
each  having  a  thickness  of  4.0  cm,  were  made  of  transparent 
Lucite  material.  Through  the  transparent  enclosure  plate, 
two-phase  flow  characteristics  of  gas  CO2  and  methanol 
solution  in  the  anode  flow  field  could  be  distinctly  visu¬ 
alized  and  recorded  by  the  image  recording  system.  It  is 
worth  pointing  out  that  the  transparent  Lucite  material, 
having  a  well-insulated  property,  can  effectively  keep  the 
heat  in  the  ME  A  and  the  bipolar  plates  from  dissipating  to 
the  surroundings  under  the  condition  of  higher  operation 
temperatures. 


2.2.  Test  loop 

The  experiments  were  carried  out  in  the  test  loop  shown 
in  Fig.  3.  Methanol  solution  was  driven  by  a  digital  HPLC 
micro-pump  (Series  III),  which  can  precisely  control  the 
liquid  flow  rate  from  0.01  to  10.0  ml  min-1  with  an  error 
of  2%  of  reading.  Before  entering  the  cell,  methanol  so¬ 
lution  was  pre-heated  to  a  desired  temperature  by  a  heater 
connected  to  a  temperature  controller.  The  mixture  of  gas 
CO2  and  unreacted  methanol  solution  was  drained  from  the 
cell  and  cooled  down  when  passing  through  a  cooling  sys¬ 
tem.  Gas  CO2  was  separated  from  the  mixture  and  released 
to  atmosphere,  while  the  unreacted  methanol  solution  was 
re-collected  into  a  chemical  liquid  tank.  Simultaneously, 
high  purity  oxygen  of  99.999%  as  oxidant  was  provided  to 
the  cathode  side  of  the  cell  without  humidification.  The  flow 
rate  of  oxygen  was  controlled  and  measured  by  a  mass  flow 
meter  (Omega  FMA-7 105E)  combined  with  a  multiple  chan¬ 
nel  indicator  (Omega  FMA-5876A),  which  can  control  the 
oxygen  flow  rate  from  0  to  500  standard  cubic  centimeter 
per  minute  (SCCM)  with  an  error  of  1%  of  full  scale.  Sim¬ 
ilar  to  methanol  solution,  oxygen  gas  was  also  pre-heated 
to  the  desired  temperature,  and  then  flowed  into  the  cathode 
side  of  the  cell.  Liquid  or  gaseous  water  would  be  produced 
depending  on  the  oxygen  flow  rate  and  cell  working  temper¬ 
ature.  With  the  help  of  another  cooling  system,  water  was 
condensed  and  collected  to  the  chemical  liquid  tank,  and  the 
excess  oxygen  gas  was  released  to  atmosphere. 

In  addition,  a  purge  system  and  a  degas  system  were  built 
in  the  present  test  loop.  Prior  to  each  new  test,  the  gas  purge 
system  was  used  to  sweep  out  the  residual  fuel  and  oxidant 


Methanol 

Solution 

Reservoir 


Produced, r  «r  Unreacted 
HjO  Methanol 

Solution 


Fig.  3.  Schematic  of  the  DMFC  test  loop. 
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Fig.  4.  Cell  orientations:  (a)  horizontal  orientation  with  the  anode  below 
the  ME  A  (i 9  =0°);  (b)  vertical  orientation  (6  =  90°);  (c)  horizontal 
orientation  with  the  anode  above  the  MEA  (i 9  =  180°). 


in  the  cell  using  nitrogen  gas.  In  order  to  prevent  methanol 
solution  from  exposing  to  atmosphere,  resulting  in  fuel  dete¬ 
rioration,  the  degas  system  was  used  to  inject  helium  gas  into 
the  closed  methanol  solution  tank  during  the  experiments. 

2.3.  Image  recording  system 

A  video  and  photo  imaging  system,  including  a  JVC  CCD, 
a  NAVITAR  6X  CCD  C-Mount  Len  and  a  digital  camcorder 
(SONY  DCR-TRV900E),  was  employed  to  capture  the  im¬ 
ages  of  the  two-phase  flow  characteristics  of  gas  CO2  and 
methanol  solution  in  the  anode  flow  field.  Being  zoomed  in 
by  the  NAVITAR  6X  CCD  C-Mount  Lens,  the  images  of  the 
anode  flow  field  were  shot  by  the  JVC  CCD  and  recoded 
on  a  mini  digital  video  cassette  by  the  digital  camcorder. 
A  shutter  speed  of  1/10,000  s  and  a  recording  speed  of  25 
frames/s  were  selected  to  visualize  and  record  each  flow  pat¬ 
tern.  A  600-W  spotlight  (ARRILITE  600)  was  employed  to 
meet  the  lighting  requirement  for  capturing  images  at  the 
high  shutter  speed  of  1/10,000  s.  The  bubble  velocity  was 
measured  from  consecutive  images  of  the  bubble  replayed  in 
slow  motion,  and  was  calculated  based  on  the  time  interval 
of  1/25  s  for  the  bubble  to  travel  a  given  distance. 

2.4.  Electrochemical  instrumentation  and  test  conditions 

In  this  work,  an  Arbin  BT2000  electro-load  interfaced  to 
a  computer  was  employed  to  control  the  cell  operation  con¬ 
dition  and  to  measure  voltage-current  (polarization)  curve. 
A  pre- written  schedule  program  was  run  to  realize  the  con¬ 
trol  and  measurement  for  a  desired  test  condition.  In  this 
schedule  program,  the  control  mode  such  as  current  or  volt¬ 
age  discharge,  the  accurate  discharging  value  and  duration, 
and  data-sampling  period  were  set  in  detail.  All  the  tests  in 
this  work  were  conducted  by  selecting  current  discharge  as 
control  mode  and  scanning  the  current  value  from  zero  (i.e., 
open  circuit)  to  a  possible  maximum  value  with  an  incre¬ 
ment  of  0.16  A  (i.e.,  10  mA  cm-2  for  the  present  MEA  area 
of  4.0  x  4.0  cm2)  each  step. 

All  the  experiments  reported  in  this  work  were  performed 
at  the  same  cathode  operation  condition,  i.e.  a  constant  oxy¬ 
gen  gas  flow  rate  of  100  SCCM  and  atmospheric  pressure 
(0.1  MPa).  As  shown  in  Fig.  4,  the  experiments  were  per¬ 
formed  for  three  different  cell  orientations:  vertical  orienta¬ 


tion  ( 0  =  90°),  horizontal  orientation  with  the  anode  below 
the  MEA  (0  =  0°),  and  horizontal  orientation  with  the  anode 
above  the  MEA  ( 0  =  180°). 


3.  Results  and  discussion 


3.1.  General  observation  of  the  CO2  gas  bubble  behavior 


The  images  of  the  CO2  gas  bubble  behavior  in  the  anode 
flow  field  corresponding  to  various  operation  conditions  are 
shown  in  Figs.  6,  8,  12  and  16.  Thanking  to  the  reflection 
of  liquid-gas  interface,  the  CO2  bubble  formation,  growth, 
and  departure  from  the  surface  of  the  gas  diffusion  layer  and 
their  motions  in  the  flow  channel  were  distinctly  captured 
by  the  image  recording  system.  As  can  be  seen  from  these 
images,  with  the  help  of  lighter  liquid-gas  interface,  gas 
bubbles  can  be  easily  distinguished  from  liquid  methanol 
solution.  The  region  enveloped  by  a  lighter-grey  liquid-gas 
interface  represents  a  gas  area,  whereas  the  rest  region  in 
the  flow  field  is  the  liquid  area. 

In  this  work,  we  observed  that  the  formation  of  CO2  gas 
bubbles  was  not  uniform  over  the  surface  of  the  gas  diffu¬ 
sion  layer:  at  some  sites,  gas  bubbles  were  generated  contin¬ 
uously,  whereas  at  the  other  sites,  no  bubble  occurred  at  all. 
This  non-uniform  distribution  of  bubbles  is  attributed  to  the 
fact  that  both  the  diffusion  layer  and  the  catalyst  layer  can¬ 
not  be  fabricated  with  perfect  uniform  surfaces.  At  each  ac¬ 
tive  nucleation  site,  we  found  that  gas  bubble  formed,  grew, 
and  departed  from  the  surface  of  the  diffusion  layer  with 
the  help  of  the  crosscurrent  flowing  methanol  solution.  As 
such,  for  a  given  operation  condition,  a  process  of  bubble 
formation  and  departure  repeated  periodically.  After  depart¬ 
ing  from  the  surface  of  the  gas  diffusion  layer,  gas  bubbles, 
usually  small  and  in  a  spherical  shape,  traveled  upwards  to 
the  upper  surface  of  the  channel  due  to  buoyancy,  and  moved 
towards  the  outlet  of  the  flow  field  with  liquid  methanol  so¬ 
lution.  Traveling  along  the  flow  channel  towards  the  outlet, 
the  bubbles  gradually  grew  up  due  to  coalescing  with  the 
newly  produced  small  bubbles  from  the  surface  of  the  gas 
diffusion  layer. 

Fig.  5  presents  the  cell  polarization  behavior  and  the 
power  output  for  the  case  when  the  cell  was  oriented  verti¬ 
cally  ( 0  =  90°)  and  with  1  M  methanol  solution  supplied  at 
a  flow  rate  of  1.0  ml  min-1  at  60  °C.  The  images  of  the  CO2 
bubble  behavior  for  selected  current  densities,  correspond¬ 
ing  to  Point  A  to  F  in  Fig.  5,  are  shown  in  Fig.  6.  A  progres¬ 
sive  scan  of  the  six  images  shown  in  Fig.  6  indicates  that  the 
quantity  of  gas  CO2  bubbles  in  the  flow  field  increased  pro¬ 
gressively  as  the  current  density  was  increased  from  10  to 
200mA cm-2.  The  increase  in  the  quantity  of  CO2  bubbles 
with  increasing  the  current  density  can  be  readily  explained 
from  the  following  equation: 


dco2  = 


AMco2  . 

- -1 

6Fpco2 
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Fig.  5.  Cell  performance  for  the  vertical  orientation  ( 6  =  90°)  with  1M  methanol  solution  fed  at  1.0  ml  min  1  and  at  60  °C. 


where  qco2  represents  the  volume  of  gas  CO2  and  A  is 
the  MEA  area.  At  low  current  densities  (e.g.  10  mA  cm-2), 
as  can  be  observed  from  the  image  for  Point  A,  a  rather 
small  number  of  discrete  gas  bubbles  were  generated  and 
present  in  the  flow  field.  Therefore,  at  low  current  densi¬ 
ties,  the  two-phase  flow  pattern  can  be  regarded  as  bubbly 
flow  [29]  over  the  entire  flow  field.  It  is  also  worth  men¬ 
tioning  that  when  the  current  density  was  sufficiently  low 
(less  than  5  mA  cm-2  for  the  cell  in  this  work),  no  CO2  gas 
bubble  was  observed  over  the  entire  flow  field.  This  phe¬ 
nomenon  is  attributed  to  the  fact  that  rather  small  quantity 


of  gas  CO2  was  produced,  which  might  completely  be  dis¬ 
solved  in  the  aqueous  methanol  solution  [30].  With  the  in¬ 
crease  in  the  current  density,  as  seen  from  the  images  for 
Point  B  and  C,  some  of  small  spherical  bubbles  grew  up 
with  a  fast  growth  rate  and  eventually  became  slug  bubbles, 
which  have  a  semi-spherically  shaped  head  and  tail,  flat  at 
top  and  bottom  boundaries.  Unlike  the  small  spherical  bub¬ 
bles,  slug  bubbles  usually  were  long  and  spanned  the  entire 
channel  cross-section.  As  a  result,  bubbly  flow  transitioned 
to  slug  flow  [29]  in  the  middle  region  of  the  flow  field,  and 
bubbly  flow  and  slug  flow  coexisted  in  the  flow  field.  With 


Fig.  6.  CO2  bubble  behavior  at  different  current  densities  for  the  vertical  orientation  (1 9  =  90°)  with  1  M  methanol  solution  fed  at  1.0  ml  min  1  and  at  60  °C. 
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the  further  increase  in  the  current  density,  as  seen  from  the 
images  for  Point  D  to  F,  slug  flow  was  present  almost  in 
entire  flow  field,  and  bubbly  flow  was  only  observed  in  the 
entrance  region  of  the  flow  field.  It  is  important  to  note  from 
the  image  for  Point  F  that  at  high  current  densities,  rather 
long  gas  slugs  formed  and  occupied  almost  the  entire  chan¬ 
nel  cross  section,  causing  the  effective  contact  area  between 
the  liquid  fuel  and  the  gas  diffusion  layer  to  be  extremely 
small.  Under  such  a  situation,  the  long  gas  slugs  may  restrict 
the  continuous  supply  of  methanol  through  the  gas  diffusion 
layer  to  the  catalyst  surface,  eventually  leading  to  deteriora¬ 
tion  of  methanol  mass  transfer.  During  the  experiments,  we 
also  found  that  for  the  cases  when  the  current  density  ex¬ 
ceeded  200  mA  cm-2,  the  CO2  bubble  behavior  in  the  flow 
field  exhibited  no  significant  difference  as  compared  with  the 
case  for  i  =  200mA  cm-2.  For  liquid-gas  two-phase  flows, 
the  variation  in  flow  patterns  is  rather  sensitive  to  the  change 
in  the  gas  void  fraction  at  intermediate  values.  However,  as 
the  gas  void  fraction  is  sufficiently  high  or  low  (high  or  low 
current  densities),  the  change  in  the  flow  patterns  becomes 
less  sensitive  to  the  gas  void  fraction.  In  addition,  at  high 
current  densities,  a  further  increase  in  the  quantity  of  gas 
CO2  only  leads  to  a  negligible  change  in  void  fraction  due 
to  the  increase  in  the  gas-liquid  slip  ratio  [31]. 

It  is  worth  mentioning  that  for  the  present  flow  field  con¬ 
sisting  of  a  single  serpentine  channel,  the  channel-blocking 
phenomenon  caused  by  CO2  gas  slugs  was  never  found  un¬ 
der  all  the  test  conditions  in  this  work,  although  this  phe¬ 
nomenon  was  found  in  the  flow  field  consisting  of  parallel 
channels  in  the  work  by  Argyropoulos  et  al.  [22].  For  the 
present  flow  field  design  with  a  single  serpentine  channel, 
there  is  no  flow  dead  end;  for  a  given  flow  rate,  gas  slugs 
may  cause  an  increase  in  the  pressure  drop  but  can  never 
block  the  channel. 


3.2.  Effect  of  cell  orientations 

The  effect  of  cell  orientations  on  the  cell  performance  for 
the  case  when  1  M  methanol  solution  was  fed  at  0.5  ml  min- 1 
and  at  60  °C  is  shown  in  Fig.  7.  It  is  clear  from  this  fig¬ 
ure  that  the  cell  orientation  has  a  significant  effect  on  the 
cell  performance,  especially  in  the  mass  transport  limitation 
region:  the  best  performance  was  achieved  when  the  cell 
was  orientated  vertically  (0  =  90°);  the  performance  became 
worse  when  the  cell  was  orientated  horizontally  with  either 
the  anode  below  the  ME  A  (6  =  0°)  or  the  anode  above  the 
ME  A  (0  =  180°).  Let  us  now  turn  our  attention  to  the  CO2 
gas  bubble  behavior  for  different  cell  orientations  at  a  se¬ 
lected  current  density  (150  mA  cm-2),  as  shown  in  Fig.  8. 
It  is  found  from  Fig.  8  that  more  long  gas  slugs  formed  in 
the  flow  field  for  the  cases  of  0  =  0°  and  0  =  180°,  whereas 
for  the  case  of  0  =  90°,  the  number  of  the  long  gas  slugs 
was  reduced.  The  effect  of  the  cell  orientation  on  the  CO2 
gas  bubble  behavior  is  attributed  to  buoyancy  force.  For  the 
vertical  orientation,  buoyancy  force  assists  the  removal  of 
CO2  gas  bubbles  on  the  top  of  the  forced  circulation  of 
methanol  solution,  thereby  tending  to  suppress  the  formu¬ 
lation  of  long  gas  slugs.  However,  as  the  cell  was  oriented 
horizontally,  the  effect  of  buoyancy  became  small,  tending 
to  accumulate  gas  bubbles  in  the  flow  field  and  eventually 
to  form  long  gas  slugs.  We  estimated  the  bubble  velocity 
by  analyzing  the  recorded  video  and  found  that  gas  bubbles 
moved  faster  for  the  case  of  9  =  90°  than  for  the  cases  of 
0=0°  and  6  =  180°.  For  instance,  at  the  current  density 
of  150  mA  cm-2,  the  bubble  velocity  in  the  exit  region  was 
about  55.7  mms-1  for  the  case  of  0  =  90°,  whereas  it  was 
reduced  to  48.7  and  50.0 mms-1  for  the  cases  of  0  =  0° 
and  6  =  180°,  respectively.  Therefore,  it  can  be  concluded 
that  the  higher  cell  performance  for  the  case  of  the  vertical 


V,  0  =  90°). 
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(a)  0=0°  (b)  0=180° 


Fig.  8.  CO2  bubble  behavior  for  different  cell  orientations  at  current  density  of  150  mA  cm  2  and  with  1M  methanol  solution  fed  at  0.5  ml  min  1  and 
at  60  °C. 


orientation  shown  in  Fig.  7  is  due  to  buoyancy  force,  which 
assists  the  removal  of  gas  bubbles  in  the  flow  field. 

It  is  also  interesting  to  note  from  Fig.  7  that  the  horizontal 
cell  orientation  with  the  anode  above  the  MEA  (0  =  180°) 
gave  better  performance  than  with  the  anode  below  the  MEA 
(0  =0°).  This  phenomenon  can  be  explained  as  follows. 
For  the  horizontal  cell  orientation  with  the  anode  above  the 
MEA  ( 6  =  180°),  bubbles  generated  from  the  lower  catalyst 
layer,  penetrated  the  diffusion  layer,  and  traveled  upwards  to 
the  upper  surface  of  the  flow  channel  (the  transparent  enclo¬ 
sure)  due  to  buoyancy.  Under  this  situation,  more  gas  bub¬ 
bles  accumulated  near  the  upper  surface  of  the  flow  channel 
but  did  not  hinder  the  access  of  methanol  solution  to  the  cat¬ 
alyst  layer.  However,  for  the  horizontal  cell  orientation  with 
the  anode  below  the  MEA  ( 0  =  0°),  generated  bubbles  will 
accumulate  near  the  upper  surface  of  the  flow  channel  (the 
diffusion  layer),  hindering  the  access  of  methanol  solution 
to  the  catalyst  layer  and  lowering  the  cell  performance. 

Figs.  9  and  10  show  the  effect  of  the  cell  orientations  on 
the  cell  performance  for  the  case  when  the  methanol  solu¬ 


tion  flow  rate  was  increased  to  1.0  and  8.0  ml  min-1  while 
the  other  parameters  were  kept  the  same  as  the  case  shown 
in  Fig.  7.  The  behavior  similar  to  Fig.  7  can  also  be  observed 
from  Figs.  9  and  10:  the  best  performance  was  achieved  for 
the  case  of  0  =  90°,  whereas  the  worst  cell  performance  was 
obtained  for  the  case  of  0  =0°.  However,  with  increasing 
the  methanol  solution  flow  rate,  the  effect  of  the  cell  orien¬ 
tations  on  the  cell  performance  became  weaker,  because  the 
forced  circulation  furnished  by  the  pump  became  larger  than 
buoyancy  at  higher  flow  rates.  This  finding  suggests  that  at 
low  methanol  flow  rates,  cautions  have  to  be  taken  to  the 
cell  orientation  to  ensure  high  cell  performance. 

3.3.  Effect  of  the  methanol  solution  flow  rate 

The  effect  of  the  methanol  solution  flow  rate  on  the 
cell  performance  at  a  low  operating  temperature  of  23  °C 
is  shown  in  Fig.  11.  The  experiments  were  conducted  for 
the  vertically-orientated  cell  (0  =  90°)  with  1  M  methanol 
solution  fed  at  flow  rates  of  0.25,  0.5,  1.0,  2.0,  4.0,  and 


Fig.  9.  Effect  of  cell  orientation  on  the  cell  performance  with  1M  methanol  solution  fed  at  1.0  ml  min  1  and  at  60  °C  (O,  6  =  0°;  □,  6  =  180°; 
V,  0  =  90°). 
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V,  0  =  90°). 


Fig.  11.  Effect  of  the  methanol  solution  flow  rate  on  the  cell  performance  for  the  vertical  orientation  (6  =  90°)  and  using  1  M  methanol  solution  at  23  °C 
(■,  0.25  ml  min-1,  O,  0.5  ml  min-1,  □,  1.0  ml  min-1,  v>  2.0  ml  min-1,  ,  4.0  ml  min-1,  ,  8.0  ml  min-1). 


8.0  ml  min-1.  As  seen  from  Fig.  11,  an  increase  in  the 
methanol  solution  flow  rate  from  0.25  to  1.0  ml  min-1  led 
to  the  improved  cell  performance.  However,  the  cell  perfor¬ 
mance  became  worse  as  the  flow  rate  was  further  increased 
from  1.0  to  8.0  ml  min-1.  This  can  be  explained  by  the  rep¬ 
resentative  images  of  the  gas  CO2  bubble  behavior  for  the 
current  density  of  60  mA  cm-2  as  shown  in  Fig.  12(a-f), 
corresponding,  respectively,  to  the  methanol  flow  rates  of 
0.25,  0.5,  1.0,  2.0,  4.0,  and  8.0  ml  min-1.  A  progressive 
scan  of  the  six  images  in  the  Fig.  12  indicates  that  both  the 
average  bubble  size  and  the  quantity  of  CO2  gas  bubbles  in 
the  flow  field  were  gradually  reduced  with  increasing  the 
methanol  solution  flow  rate  from  0.25  to  8.0  ml  min-1.  For 


lower  methanol  solution  flow  rates  of  0.25  and  0.5  ml  min-i 
shown  in  Fig.  12a  and  b,  some  long  gas  slugs  were  present 
in  the  flow  field,  which  reduced  the  effective  contact  area 
between  methanol  and  the  gas  diffusion  layer  and,  thus  led 
to  poor  cell  performance.  As  the  methanol  solution  flow 
rate  was  increased  to  1.0  ml  min-1,  as  shown  in  Fig.  12c, 
the  number  of  gas  CO2  bubble  was  reduced  and  the  long 
gas  slugs  disappeared  from  the  flow  field  because  of  the  fact 
that  the  sweeping  rate  of  the  gas  bubbles  was  increased  with 
the  methanol  solution  flow  rate.  As  a  result,  an  increased 
effective  contact  area  between  liquid  methanol  and  the  gas 
diffusion  layer  enhanced  the  mass  transfer  of  methanol  and 
hence,  improved  the  cell  performance  as  seen  in  Fig.  11. 
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Fig.  12.  CO2  bubble  behavior  for  different  flow  rates  at  60mA cm  2,  6  =  90°  and  using  1M  methanol  solution  at  23  °C. 


For  the  higher  methanol  solution  flow  rates  of  2.0,  4.0,  and 
8.0  ml  min-1,  as  shown  in  Fig.  12(d-f),  the  quantity  of  gas 
CO2  bubble  decreased  with  increasing  methanol  solution 
flow  rates.  Under  such  a  situation,  the  mass  transfer  of 
methanol  was  not  a  serious  problem  any  more.  However, 
an  increase  in  the  methanol  solution  flow  rate  is  accompa¬ 
nied  by  an  increase  in  the  static  pressure  in  the  flow  field. 
A  higher  static  pressure  tends  to  increase  in  the  methanol 
crossover  from  the  anode  to  the  cathode,  leading  to  deterio¬ 
ration  in  the  cell  performance.  As  seen  from  Fig.  11,  the  cell 


performance  became  worse  and  worse  when  the  methanol 
solution  flow  rate  was  increased  from  1.0  to  8.0  ml  min-1. 

Figs.  13  and  14  show  the  effect  of  the  methanol  solu¬ 
tion  flow  rates  of  0.5,  1.0,  2.0,  4.0,  and  8.0  ml  min-1  on  the 
cell  performance  when  temperature  was  increased  to  40  and 
60  °C,  but  the  other  parameters  were  kept  the  same  as  the 
case  shown  in  Fig.  11.  It  is  seen  that  at  a  moderate  oper¬ 
ating  temperature  of  40  °C,  the  best  cell  performance  was 
achieved  at  a  flow  rate  of  1.0  ml  min-1,  but  either  a  lower 
flow  rate  (0.5  ml  min-1)  or  a  higher  flow  rate  (2.0,  4.0,  and 


Fig.  13.  Effect  of  the  methanol  solution  flow  rate  on  the  cell  performance  for  6  =  90°  and  using  1M  methanol  solution  at  40  °C  (O,  0.5  ml  min  1 ,  □, 
1.0  ml  min-1,  v>  2.0mlmin-1,  ,  4.0mlmin-1,  ,  8.0mlmin-1). 
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Fig.  14.  Effect  of  the  methanol  solution  flow  rate  on  the  cell  performance  for  6  =  90°  and  using  1M  methanol  solution  at  60  °C  (O,  0.5  ml  min  ,  □, 
1.0  ml  min-1,  v>  2.0mlmin-1,  ,  4.0mlmin-1,  ,  8.0mlmin-1). 


8.0  ml  min-1)  led  to  a  deterioration  in  the  cell  performance. 
Similarly,  at  a  high  operating  temperature  of  60  °C,  an  in¬ 
crease  in  flow  rate  from  0.5  to  2.0  ml  min-1  improved  the 
cell  performance  significantly,  but  a  further  increase  from 
2.0  to  8.0  ml  min-1  resulted  in  a  deterioration  in  the  cell  per¬ 
formance.  Furthermore,  a  progressive  scan  of  the  cell  perfor¬ 
mance  curves  shown  in  Figs.  11,  13  and  14,  corresponding, 
respectively,  to  temperature  of  23,  40,  and  60  °C,  indicates 
that  the  effect  of  the  methanol  solution  flow  rate  on  the  cell 
performance  became  less  sensitive  with  increasing  the  cell 
operating  temperature.  These  experimental  results  revealed 
that  the  methanol  solution  flow  rate  has  a  significant  effect 
on  the  cell  performance,  especially  at  low  operating  temper¬ 
atures.  This  finding  is  significant  for  portable  DMFCs  that 


operate  at  room  temperature.  Cautions  have  to  be  taken  to 
ensure  an  optimal  flow  rate  to  be  designed  for  those  DMFCs 
operating  at  low  temperatures. 

3.4.  Effect  of  cell  operating  temperature 

Fig.  15  shows  the  effect  of  temperatures  on  the  cell  per¬ 
formance  for  0  =  90°  and  with  1  M  methanol  solution  fed  at 
2.0  ml  min-1 .  As  seen  from  this  figure,  the  cell  performance, 
including  OCV,  potential  and  power  density  increased  with 
increasing  temperature  from  23  to  60  °C.  The  maximum 
power  density  was  about  55  mW  cm-2  at  60  °C.  Fig.  16 
shows  the  images  of  gas  CO2  bubble  behavior  at  a  selected 
current  density  of  80  mA  cm-2  corresponding  to  the  case 


Fig.  15.  Effect  of  temperature  on  the  cell  performance  for  6  =  90°  and  with  1M  methanol  solution  fed  at  2.0  ml  min  1  (O,  23  °C,  □,  40  °C,  v»  60  °C). 
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Fig.  16.  Effect  of  temperature  on  the  CO2  bubble  behavior  for  0  =  90°  at  80  mA  cm 


with  1M  methanol  solution  fed  at  2.0  ml  min 


in  the  Fig.  15.  As  can  be  seen  from  the  images  shown  in 
Fig.  16(a-c),  the  quantity  of  CO2  gas  bubbles  in  the  flow 
field  increased  gradually  with  increasing  the  cell  operating 
temperature  from  23  to  60  °C.  As  indicated  by  Eq.  (4),  the 
quantity  of  gas  CO2  depends  on  the  current  density  for  a 
given  temperature.  However,  the  images  shown  in  Fig.  16 
indicate  that  the  quantity  also  depends  on  the  cell  tempera¬ 
ture  at  a  fixed  current  density.  There  are  two  reasons  associ¬ 
ated  with  this  phenomenon.  First  for  a  given  static  pressure 
(depending  on  the  methanol  solution  flow  rate),  the  CO2 
solubility  in  the  aqueous  methanol  solution  decreases  with 
increasing  temperature,  thereby  resulting  in  more  CO2  gas 
bubbles  in  the  flow  field  at  a  higher  temperature.  Secondly, 
based  on  the  ideal  gas  equation  of  state,  the  specific  volume 
of  CO2  increases  with  temperature,  also  leading  to  the  pres¬ 
ence  of  a  larger  volume  of  gas  CO2  in  the  flow  field.  As  a 
matter  of  fact,  with  the  aid  of  the  ideal  gas  equation  of  state: 

Pco2  =  pMCo2/RT  (5) 


Eq.  (4)  can  be  recast  as: 


q  co2  = 


AR 

- iT 

6Fp 


Based  on  Eq.  (6),  it  can  be  estimated  that  for  the  ME  A 
area  of  4.0  x  4.0  cm2  in  this  work  at  a  given  current  density, 
the  increase  in  temperature  from  23  to  60  °C  would  cause 
an  increase  in  the  volume  of  gas  CO2  by  about  12.5%. 

It  is  also  interesting  to  note  from  Fig.  16  that  for  a  fixed 
current  density,  the  average  size  of  the  bubbles  in  the  flow 
field  increased  with  temperature.  This  phenomenon  is  related 
to  surface  tension  of  methanol  solution,  which  decreases 
with  temperature.  Thus,  small  bubbles  tended  to  coalesce 
with  each  other  and  to  form  larger  bubbles.  In  addition,  the 
evaporation  of  methanol  solution  at  a  higher  temperature 
might  also  contribute  to  this  phenomenon.  As  can  be  seen 
in  Fig.  16(a-c)  at  the  current  density  of  80  mA  cm-2,  gas 
slugs  gradually  became  predominated  in  the  flow  field  with 
increasing  temperature.  When  the  temperature  reached  to 
60  °C,  plenty  of  gas  slugs  formed  and  led  to  a  reduced  effec¬ 
tive  contact  area  between  liquid  methanol  solution  and  the 
gas  diffusion  layer,  causing  a  negative  effect  on  the  cell  per¬ 
formance.  However,  an  increase  in  temperature  not  only  in¬ 
creases  the  activity  of  the  catalysts,  but  also  leads  to  a  higher 


diffusion  coefficient  of  methanol,  both  of  which  are  favor¬ 
able  for  improving  the  cell  performance.  This  indicates  that 
the  increase  in  both  the  catalyst  activity  and  mass  transfer 
coefficient  of  methanol  caused  by  the  increase  in  tempera¬ 
ture  suppress  the  negative  effect  of  large  gas  slugs  formed 
at  high  temperatures. 


4.  Conclusions 

A  visual  investigation  of  CO2  gas  bubble  behavior  in¬ 
side  the  anode  flow  field  consisting  of  a  single  serpentine 
channel  has  been  performed.  The  experiments  have  shown 
that  in  the  anode  flow  field  there  exist  various  two-phase 
flow  patterns  of  gas  CO2  and  methanol  solution  that  vary 
with  the  cell  operating  conditions.  The  quantity  of  gas  CO2 
bubbles  in  the  flow  field  increased  with  the  current  density. 
At  low  current  densities,  small  discrete  bubbles  appeared  in 
the  flow  field.  At  moderate  current  densities,  a  number  of 
gas  slugs  formed  in  addition  to  small  discrete  bubbles.  And 
at  high  current  densities,  the  flow  field  was  predominated 
by  rather  long  gas  slugs.  It  has  been  shown  that  the  pres¬ 
ence  of  long  gas  slugs  in  the  flow  field  led  to  a  determina¬ 
tion  in  the  cell  performance.  The  experiments  also  indicate 
that  the  cell  orientation  has  a  significant  effect  on  the  cell 
performance,  especially  at  low  methanol  flow  rates;  for  the 
present  flow  field  design  the  best  cell  performance  can  be 
obtained  when  the  cell  is  orientated  vertically.  This  visually 
study  shows  that  an  increase  in  the  methanol  solution  flow 
rate  might  tend  to  reduce  the  average  size  and  the  number  of 
gas  slugs  in  the  flow  field,  thereby  enhancing  the  methanol 
mass  transfer.  However,  high  flow  rates  may  lead  to  an  in¬ 
creased  methanol  crossover,  resulting  in  deterioration  in  the 
cell  performance.  It  has  also  been  revealed  that  the  effect  of 
the  methanol  solution  flow  rate  became  more  pronounced  at 
low  temperatures.  Therefore,  for  portable  DMFCs  that  op¬ 
erate  at  room  temperature,  special  caution  has  to  be  taken 
for  determining  an  optimal  flow  rate.  The  study  on  the  ef¬ 
fect  of  temperatures  reveals  that  the  cell  performance  got 
improved  with  increasing  temperature,  although  large  gas 
slugs  formed  in  the  flow  field  at  high  temperature.  Further¬ 
more,  for  the  present  flow  field  consisting  of  a  single  ser¬ 
pentine  channel,  channel-blocking  phenomenon  caused  by 
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CO2  gas  slugs  were  never  found  for  all  the  test  conditions 

in  this  work. 
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